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1. Introduction 

The model of a random walk (RW) provides a good description of diffusion processes, 
such as for example encountered for electrons in metals or colloidal particles in solution 
[1]. The averaged mean square displacement of the diffusive particle at time t (or, 
equivalently, after t steps on a lattice) scales as 

where in a non-disordered medium i/rw = 1/2, independently of the space dimension 
d. A RW is a fractal object, with fractal dimension c^rw = ^/^rw- The number of all 
possible trajectories Ct for a randomly walking particle of t steps can be found exactly: 
Ct = Zq, where zq is the coordination number of the corresponding lattice. 

Forbidding the trajectory of a random walk to cross itself, we obtain a self- 
avoiding walk (SAW), which is one of the most successful in describing the universal 
configurational properties of a long, flexible single polymer chain in good solvent [2]. 
The average squared end-to-end distance (i?^) and the number of configurations Cn of 
SAWs with A^ steps on the underlying lattice obey the scaling laws: 

(R^) ~ Ar2^sAw^ Cn ~ ^^A^tsaw-i^ ^2) 

where t'sAW) 7saw are universal exponents that only depend on the space dimensionality 
d, and z is a non-universal fugacity, counting the average number of accessible nearest- 
neighbor sites. The properties of SAWs on a regular lattice have been studied in detail 
both in analytical approaches [3l HI O [6] and computer simulations [3 [HJ [9l [TOl [HI [12] . 
For example, in the space dimension d=3 one finds within the frame of the field- 
theoretical renormalization group approach i/sAw=0.5882 ± 0.0011 [6j and Monte Carlo 
simulations give z/saw=0.5877 ± 0.0006 [H]. For space dimensions d above the upper 



critical dimension (iup=4, the effect of self-avoidance becomes irrelevant and SAWs 
behave effectively as random walks with exponents i/rw = 1/2, 7rw = 1- 

The problem of random walks in disordered media is of great interest since it is 
connected with a large amount of physical phenomena: transport properties in fractures 
and porous rocks, the anomalous density of states in randomly diluted magnetic systems, 
silica aerogels and in glassy ionic systems, diffusion-controlled fusion of excitations in 
porous membrane films etc. (see, e.g., Ref. [13j for a review). Similarly, SAWs on 
randomly diluted lattices may serve as a model of linear polymers in a porous medium. 



Table 1. Critical concentration pc of site-diluted lattices and fractal dimensions of 
percolation clusters, df., and the backbone of percolation clusters, dp^, for different 
space dimensions d. 

F jB 



Pc dp, d. 



2 0.592746 [Hj 91/49 [Hj 1.650 ± 0.005 [ISj 

3 0.31160 dH 2.51 ± 0.02 |T7] 1.86 ± 0.01 [H] 

4 0.19688 JB| 3.05 ±0.05 HI] 1.95 ± 0.05 [H] 
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Much of our understanding of disordered systems comes from percolation theory 
|19j . A disordered medium can be modelled as randomly diluted lattice, with a 
given concentration p of lattice sites allowed for walking. Most interesting is the 
case, when p equals the critical concentration pc, the site-percolation threshold (see 
Table [1]), and an incipient percolation cluster can be found in the system. Studying 
properties of percolative lattices, one encounters two possible statistical averages. In 
the first, one considers only percolation clusters with linear size much larger than 
the typical length of the physical phenomenon under discussion. The other statistical 
ensemble includes all the clusters which can be found in a percolative lattice. For 
the latter ensemble of all clusters, the walks can start on any of the clusters, and 
for an A^-step walk, performed on the ith cluster, we have {R"^) ~ If, where /j 
is the averaged size of the ith cluster. In what follows, we will be interested in 
the former case, when trajectories of walks reside only on the percolation cluster. 
In this regime, the scaling laws ([1]), ([2]) hold with new exponents u^ ^ v-^^ 

[asiEiiEHiiMiiiniiiiiiiiiisiiiaiiaiiniiiTiiiHlia a 

hint to the physical understanding of this phenomenon is given by the fact that weak 
disorder does not change the dimension of a lattice, whereas the percolation cluster 
itself is a fractal object with fractal dimension rff^ dependent on d (see Table [1]). In 
this way, scaling law exponents of residing walks change with the dimension d^^ of 
the (fractal) lattice on which the walk resides. Since (iup=6 for percolation [19], the 
exponents z/gAw(^ > 6)=l/2, -il\y^[d > 6)=1. 

Our present paper aims to supplement the studies of random and self-avoiding walks 
on percolative lattices by obtaining numerical values for exponents, governing the scaling 
behavior of the models, up to rf = 4 by computer simulations. The layout of the paper is 
as follows: in the next section, we will present in detail the procedure of extracting the 
percolation cluster and its backbone on disordered lattices at the percolation threshold. 
In section III we describe the pruned-enriched Rosenbluth algorithm, applied to study 
the scaling of self-avoiding walks, and present the results obtained. In the next section 
we consider the method for studying random walks on percolation clusters. In Section 
V, we end up by giving conclusions and an outlook. 

2. Construction of percolation cluster 

We consider site percolation on a regular lattice of edge length L = 400, 200, 50 in 
dimensions d = 2,3,4, respectively. Each site of the lattice is assigned to be occupied 
with probability pc (values of critical concentration in different dimensions are given 
in the Table [1]), and empty otherwise. To describe the procedure of extracting the 
percolation cluster, let us consider schematically the two-dimensional case. We apply 
an algorithm based on the one proposed by Hoshen and Kopelman p9]. As a first step, 
a label is prescribed to each of the occupied sites. Such a labeling process is regulated, 
we start, for example, from the first "column" of the lattice, label the occupied sites 
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upwards, and then turn to the next "column", as shown in Fig. [H left. Next, we start 
the procedure of burning the occupied sites. Again, in the same order, starting from 
the bottom of the first "column" of the lattice, for each of the labeled sites (say, i), we 
check whether its nearest neighbors are also occupied or not. If yes, two possibilities 
appear: 1) The label of the neighbor is larger than the label of site i. In this case, we 
change the label of the neighbor to that of site i. 2) The label of the neighbor is smaller 
than that of i. Then, we change the label of site i to that of the neighbor. 

Such a burning procedure is applied until no more change of site labels is needed. 
As a result, we end up with groups of clusters of occupied sites with the same labels 
(Fig. [U right). Finally, we check, whether there exists a cluster, that wraps around 
the lattice in all d directions. If yes, we have found the percolation cluster (Fig. [2]). If 
not, this disordered lattice is rejected and a new one is constructed. Note, that on finite 
lattices the definition of spanning clusters is not unique (e.g., one could consider clusters 
connecting only two opposite borders), but all definitions are characterized by the same 
fractal dimension and are thus equally legitimate. The here employed definition has the 
advantage of yielding the most isotropic clusters. Note also that directly a.t p = Pc more 
than one spanning cluster can be found in the system, and the probability P{k) for at 
least k separated clusters grows with the space dimension as P{k) ~ exp{—ak'^^^'^~^^) 
[60] . In our study, we take into account only one cluster per each disordered lattice 
constructed, in order to avoid presumable correlations of the data. 

Aiming to investigate the scaling of SAWs on percolative lattices, we are interested 
in the backbone of percolation clusters, which can be defined as follows. Assume that 
each bond (or site) of the cluster is a resistor and that an external potential drop is 
applied at two ends of the cluster. The backbone is the subset of the cluster consisting 
of all bonds (or sites) through which the current flows; i.e., it is the structure left when 
all "dangling ends" are eliminated from the cluster. The SAWs can be trapped in 
"dangling ends" , therefore infinitely long chains can only exist on the backbone of the 
cluster. The algorithm for extracting the backbone of obtained percolation clusters was 
first introduced in Ref. [61] and improved in Ref. [62] • We choose the starting point 
- "seed" - at the center of the cluster, and find the chemical distance / of all the sites 
belonging to the percolation cluster to this starting point. In Fig. [3l the starting point 
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Figure 1. Procedure of site labeling and extracting the percolation cluster. 
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Figure 2. Percolation cluster on a d = 2-dimensional regular lattice of edge length 



is numbered with 0, and the chemical distance of all the other sites are depicted. The 
burning algorithm is divided into two parts. First, we start from some site of the edge 
of the lattice belonging to the percolation cluster and consider it as burning site. At the 
next step, if the nearest neighbor of this site has the chemical distance smaller than the 
burning site itself, the nearest neighbor site is burnt. Such a procedure ends when the 
"seed" site is reached. All the thus obtained burnt sites are located along the shortest 
path between the "seed" and the given site at the edge of the percolation cluster, as is 
shown in Fig. [3l The same algorithm is applied successively to all the edge sites. As a 
result, we obtain the so-called skeleton or elastic backbone |63], shown in Fig. HI left. 
In the second part of the algorithm, we consider successively each site of the elastic 
backbone, and check, whether a "loop" starts from this site. "Loop" is a path of sites, 
belonging to the percolation cluster, which is connected with the elastic backbone by no 
less than two sites. All sites of the elastic backbone together with the sites of "loops" 
form finally the geometric backbone of the cluster (see Fig. |U right). 

The results for fractal dimensions of the percolation cluster and its geometrical 
backbone in d=2, 3, 4 are compiled in Table [TJ 
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Figure 3. For all sites of a percolation cluster the chemical distance from the 
starting site is calculated. The minimal paths from all the sites on the edge of the 
percolation cluster to this starting point are found, which form the elastic backbone 
of the percolation cluster. 
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Figure 4. Elastic and geometrical backbones of the percolation cluster depicted in 
Fig.H 



3. Self-avoiding walks on percolation clusters 



3.1. The method 

For the sampling of SAWs, we use the pruned-enriched Rosenbluth method (PERM), 
proposed in the work of Grassberger [M]- The algorithm is based on ideas from the very 
first days of Monte Carlo simulations, the Rosenbluth- Rosenbluth (RR) method [7] and 
enrichment [65j- Let us consider the growing polymer chain, i.e. the nth monomer is 
placed at a randomly chosen neighbor site of the last placed (n — l)th monomer (n < N, 
where A^ is total length of polymer). The growth is stopped, if the total length A^ of the 
chain is reached, then the next chain is started to grow from the same starting point. 
In order to obtain correct statistics, if this new site is occupied, any attempt to place 
a monomer at it results in discarding the entire chain. This leads to an exponentional 
"attrition" , the number of discarded chains grows exponentially with the chain length, 
which makes the method useless for long chains. In the RR method, occupied neighbors 
are avoided without discarding the chain, but the bias is corrected by means of giving a 
weight Wn ~ (nr=i "^') ^° ^^^^ sample configuration at the nth step , where mi is the 
number of free lattice sites to place the /th monomer. When the chain of total length 
A^ is constructed, the new one starts from the same starting point, until the desired 
number of chain configurations are obtained. The configurational averaging, e.g., for 
the end-to-end distance r{N) = ^JW\N)^ then has the form: 



{t(N)) 



conf 



r'(r-w-ro)^ 



^N 



conf 



w 



N 



(3) 



where tq is the position of the starting point of the growing chains, fk is the position of 
the /cth monomer, and Z^q is the partition sum. 

The Rosenbluth method, however, also suffers from attrition: if all next neighbors 
at some step (n < A^) are occupied, i.e., the chain is running into a "dead end", 
the complete chain has to be discarded and the growth process has to be restarted. 
Combining the chain growth algorithm with population control, such as PERM (pruned- 
enriched Rosenbluth method) |64j leads to a considerable improvement of the efficiency 
by increasing the number of successfully generated chains. The weight ffuctuations of 
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the growing chain are suppressed in PERM by pruning configurations with too small 
weights, and by enriching the sample with copies of high- weight configurations |64j . 
These copies are made while the chain is growing, and continue to grow independently 
of each other. Pruning and enrichment are performed by choosing thresholds W^ and 
W^ depending on the estimate of the partition sum for the n-monomer chain. These 
thresholds are continuously updated as the simulation progresses. The zeroes iteration 
is a pure chain-growth algorithm without reweighting. After the first chain of full length 
has been obtained, we switch to W^, W^ . If the current weight Wn of an n-monomer 
chain is less than W^, a random number r = 0, 1 is chosen; if r = 0, the chain is 
discarded, otherwise it is kept and its weight is doubled. Thus, low-weight chains are 
pruned with probability 1/2. If Wn exceeds W^ , the configuration is doubled and the 
weight of each identical copy is taken as half the original weight. For a value of the 
weight lying between the thresholds, the chain is simply continued without enriching 
or pruning the sample. For updating the threshold values we apply similar rules as in 
[661 EZ]: W> = C{Zn/Zi){cn/cif and W< = 0.2W>, where c„ denotes the number of 
created chains having length n, and the parameter C controls the pruning-enrichment 
statistics. After a certain number of chains of total length A^ is produced, the given tour 
is finished and a new one starts. We adjust the pruning-enrichment control parameter 
such that on average 10 chains of total length A^ are generated per each tour [67J. Also, 
what is even more important for efficiency, in almost all iterations at least one such a 
chain was created. The number of different trajectories of SAWs with A^ steps can be 
then estimated as averaged weight: 

Cn = ]Y1 ^^"'' (4) 

conf 

where t is the number of successful tours. PERM has been applied to a wide class of 
problems, in particular study of 0-transition in homopolymers [61] , trapping of random 
walkers on absorbing lattices [68], study of protein folding [69l [67] etc. 

3.2. Results 

To study scaling properties of SAWs on percolating lattices, we have to perform two 
types of averaging: the first average is performed over all SAW configurations on a single 
percolation cluster according to ([3]); the second average is carried out over different 
realizations of disorder, i.e. over all percolation clusters constructed: 

1=1 

i=l 

where M is the number of different clusters and the index i means that a given quantity 
is calculated on the cluster i. 
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Figure 5. Disorder averaged end-to-end distance vs number of steps in double 
logarithmic scale for SAWs on the backbone of percolation clusters in d=2 (pluses), 
d=3 (stars), d=4 (squares). Lines represent linear fitting, statistical error bars are of 
the size of symbols. 



The SAW statistics crucially depends on the type of disorder averaging, namely, 
whether the disorder is considered to be "annealed" (positions of defects are in 
thermodynamical equilibrium with the system) or "quenched" (positions of defects are 



Table 2. The exponent J^g^w ^"-"^ ^ SAW on a percolation cluster. FL: 

Flory-like theories, EE: exact enumerations, RS, RG: real-space and field-theoretic 
renormalization group, MC: Monte Carlo simulations. For SAWs on the regular lattice 
one has: i^sAw(rf=2)=3/4 0], i/sAw(rf=3)=0.5882(ll) [6], z/sAw(rf > 4)=l/2. 



'^S AW \ d 


2 


3 


4 


FL Eq. Q 


7/9 


0.665 


0.594 


49J 


0.778 


0.662 


0.593 


m 


0.77 


0.66 


0.62 


m 


0.770 


0.656 


0.57 


m 


0.76 


0.65 


0.58 


EE[43] 


0.745(10) 


0.635(10) 




m 


0.770(5) 


0.660(5) 




|46| 


0.778(15) 


0.66(1) 




RS |49| 


0.778 


0.724 




ES] 


0.767 






RG P] 


0.785 


0.678 


0.595 


m 


0.796 


0.669 


0.587 


MC [3S] 


0.77(1) 






m 


0.783(3) 






141j 




0.62-0.63 


0.56-0.57 


m 


0.787(10) 


0.662(6) 




our results 


0.782 ±0.003 


0.667 ±0.003 


0.586 ±0.003 
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out of equilibrium) . As it was pointed out in Ref. [TD], the correctness of results, obtained 
in the picture of "quenched" disorder, depends on whether the location of the starting 
point of a SAW is fixed while the configurational averaging is performed, or not. In the 
latter case, one has to average over all locations and effectively this corresponds to the 
case of annealed disorder. 

An interesting question arises: what is the difference in statistics between 
SAWs walking on percolation clusters and the backbone of percolation clusters, after 
eliminating all the "dead ends"? First Kim [3H] claimed, based on a scaling argument, 
that the critical behavior on the percolation cluster is the same as that on the backbone. 

This equality was also assumed by Rammal [50] in deriving a 



z/< 



B 

SAW- 



namely i/g^w 

Flory formula for SAWs on fractal substrates. This can be easily explained: since the 
walks, which visit the dead ends are eventually terminated after a certain number of 
steps, the walks that survive in the limit N ^ oo are those confined within the backbone. 
However, in a numerical study [39] it was found, that t'^^-yv > '^saw ^^'^' moreover, that 
^SAw almost equals the value for SAWs on pure lattices. It was argued, that the averaged 
end-to-end distance of SAWs on the backbone is significantly enhanced in comparison 
with averaging on the full percolation cluster. We have checked this, comparing results 
obtained by us for the averaged end-to-end distance {r{N)) on percolation clusters and 
the backbone of percolation clusters. We conclude, that there is practically no difference 
in scaling for SAWs on both types of clusters, the SAW statistics is determined by the 
backbone of percolation clusters. 

To study the scaling properties of SAWs on the backbone of percolation clusters, we 
choose as the starting point the "seed" of the cluster, and apply the PERM algorithm, 
taking into account, that a SAW can have its steps only on the sites belonging to the 
backbone of the percolation cluster. We use lattices of size up to Linax=400, 200, 50 in 
(i=2,3,4, respectively, and performed averages over 1000 percolation clusters in each 
case. Estimates for the critical exponents t'gAw "were obtained by linear least-square 
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Figure 6. Disorder averaged distribution function rP{r, N) vs the scaling variable 
r/iV^sAw in d=2 dimensions. Lattice size L=400, number of SAW steps 7V=140 
(squares), A^=120 (triangles), A^=100 (crosses). 
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Figure 7. Disorder averaged distribution function rP{r, N) vs the scaling variable 
r/N"sAw in d—3 dimensions. Lattice size L—200, number of SAW steps iV=80 
(squares), iV=60 (triangles), iV=50 (crosses). 




Figure 8. Disorder averaged distribution function rP{r, N) vs the scaling variable 
r/iV^sAw in d=4 dimensions. Lattice size L=50, number of SAW steps A^=30 (squares), 
N=20 (triangles), 7V=15 (crosses). 



fitting (see Tables [6l [71 [8] in the Appendix). Note, that since we can construct lattices 
only of a finite size L, it is not possible to perform very long SAWs on it. For each 
L, the scaling for {r{N)) holds only up to some "marginal" number of SAWs steps 
A'marg ~ L ''^sAw j58j. Wc take this into account when analyzing the data obtained; for 
each lattice size we are interested only in values of A^ < A^arg, thus avoiding distortions, 
caused by finite-size effects. Our results for the scaling exponent t'sAw ^^^ SAWs on the 
backbone of percolation clusters [58] are given in Table [2] and compared with previous 
estimates obtained by a variety of different techniques. We see that the value of t'gAw i^ 
larger than the corresponding exponent on the pure lattice; presence of disorder leads 
to stretching of the trajectory of self-avoiding walks. 

A simple modified Flory formula for the exponent of a SAW on a percolation cluster, 
proposed a long time ago by Kremer [36] . 

(7) 



'SAW 



3/«+2) 
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Table 3. The connectivity constant z^" for a SAW on a percolation cluster. SS: 
series studies, EE: exact enumerations, MC: Monte Carlo simulations. For SAWs on 
the regular lattice one has z{d = 2) = 2.6385±0.0001 [75], z{d = 3) = 4.68404±0.00009 
[TU], zid = A) = 6.77507 ±0.00001 9 . 



zP'\d 


2 


3 


4 


SS t48j 


1.31 ±0.03 






EE[13] 


1.53 ±0.05 






m 


1.565 ±0.005 


1.462 ±0.005 




MC [39] 


1.459 ±0.003 






m 


1.456 ±0.005 


1.317 ±0.005 




our results 


1.566 ±0.005 


1.459 ±0.005 


1.340 ±0.005 


Z-Pc 


1.564 


1.460 


1.333 



gives numbers in an astonishingly good agreement with our numerical data (see Table [2]). 
For the estimates we have used the values of the fractal dimension of percolation clusters 
from Table [H Since dup = 6 for percolation and dp^{d > 6) = 4 [19], we receive from 
Eq. (jTj) that the exponent i^^Awi^ > 6) = 1/2. Note, that there exists a whole family of 
more sophisticated Flory-like theories [381 liOllSO l [STl |52l l53l[53] . 
The disorder averaged distribution function, defined via 

^=^rP(r,Ar) (8) 

r 

can be written in terms of the scaled variables r/{r) as 



rP{r, N) ~ f{r/{r)) ~ /(r/A^^'sAw). (g) 



The distribution function is normalized according to ^^, P(r, A^)=l. The numerical 
results for the distribution function in d = 2,3, and 4 are shown in Figs. [6], U\ and M 
for different N. When plotted against the scaling variable r/N^skw^ the data are indeed 
found to nicely collapse onto a single curve, using our values for the exponent t'gAw 
reported in Table [21 



Table 4. The exponent 7gAw '^'-"^ ^ SAW on a percolation cluster. FL: Flory-like 
theories, EE: exact enumerations, MC: Monte Carlo simulations. For a SAW on the 
regular lattice one has 7sAw('i = 2) = 43/32 [H], 7sAw(d = 3) = 1.1596 ± 0.0002 [5], 
7SAw((i>4) = 1. 



7sAW \ "■ 
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FL[53] 


1.384 


1.379 


1.27 


EE [43] 


1.3 ±0.1 
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1.34 ±0.05 


1.29 ±0.05 




MC [37] 


1.31 ±0.03 


1.40 ±0.02 




m 


1.26 ±0.05 


1.19 ±0.05 




our results 


1.350 ±0.008 


1.269 ±0.008 


1.250 ±0.008 
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Figure 9. Averaged connectivity constant for SAWs on the backbone of percolation 
clusters in d—2 (triangles), d=3 (squares), d—i (stars). 



Let us now turn our attention to estimating the number of different possible SAW 
configurations C^, defined by Eq. ([2]). First, let us note, that the fugacity or connectivity 
constant z is obviously affected by introducing disorder into the lattice. For the case, 
when the SAW is not confined only to the percolation cluster, namely when averaging 
over all the clusters is performed, then zP''=pcZ exactly. In Table [3] we present results of 
this estimate, taking values for p^ from Table [TJ However, since each existing bond on 
the infinite percolation cluster has a non-trivial (correlated) probability of occurrence, 
a similar argument cannot be applied to the SAWs confined to the infinite percolation 
cluster only. However, enumeration estimates [71] suggested z^^^PcZ up to Pc for SAWs 
on the percolation cluster. It turns out, that this difference from linear dependence 
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Figure 10. Disorder averaged number of SAWs configurations vs number of steps for 
SAWs on the backbone of percolation clusters in d=2 (triangles), d=3 (squares), d=A 
(stars). 
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for incipient cluster is subtle and could hardly be detected. We have estimated Zp as 
the averaged number of possibilities to perform the next step in the PERM procedure 
for SAWs on the backbone of percolation clusters (see Fig. [9]); results are presented in 
Table [31 

In the analytical study of Ref. [73], it was argued that the exponent 7, governing 
the scaling of the number of SAW configurations, is not changed by the presence of 
disorder even at p = Pc- This was supported by an exact enumeration study [i3] . 
However, this argument disagrees with results of studies [531 133 ETj, where averaging 
over single percolation clusters was performed and different values for 7gAw were found. 
In Ref. [53] it was proven, using scaling arguments, that at p = pc the exponents 7saw 
will crossover to 7gAw — 7saw + '^('^saw ~ '^saw) at p = p^, the estimates based on this 
equality are given in the first row of Table HI 

We obtained numerical estimates for 7saw' studying the behavior of the quantity 

\nC^/N = \n{A)/N + ln{z^^) + (7^^^ - 1)^^, (10) 

where A is a constant. Figure ^U\ shows these values for the backbone of percolation 
clusters in d = 2,3,4. Estimates for 7gAw ^^^ obtained by linear least-square fits (see 
Tables [9l [TOl [TTl in the Appendix). Our final results are presented in Table HI 

4. Random walks on percolation clusters 



To simulate the diffusion process in a disordered medium, the picture of the "ant in 
the labyrinth", proposed by de Gennes [E] is traditionally used. Here the walker (an 
"ant") starts at an arbitrary point on the diluted lattice and tries to move randomly 
to the nearest site. If the randomly chosen direction leads to an empty site, it moves 
and the steps increment by 1. If the chosen site is occupied by a defect (in our case. 




Figure 11. Averaged end-to-end distance vs number of steps in a double logarithmic 
scale of RWs on percolation clusters (crosses) and the backbone of percolation clusters 
(pluses) in d — 2. 
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does not belong to the percolation cluster) it stays at the current position for this time 
step. The growth is stopped, if the total number of steps A^ is performed, than the next 
trajectory is started to grow. 

After averaging the end-to-end distance over RW configurations on a single 
percolation cluster, the disorder average is carried out as in Eq. ([5]) over all constructed 
percolation clusters. Let us note that, in contrast to the SAW problem, discussed 
previously, the scaling behavior of RWs on a percolation cluster is different from that 
on its backbone. Let us remind, that statistics of long SAWs on percolation clusters is 
nevertheless determined by its backbone, since the walks, which visit the "dead ends" 
are eventually terminated after a certain number of steps. Simple random walks cannot 
be trapped in "dead ends" , however, visiting these parts of a cluster will lead to some 
"slowing down" of the diffusion process in comparison with the behavior on the backbone 
where all the dead ends are removed. This is really confirmed by analyzing our results 
for the averaged end-to-end distance of random walks on a percolation cluster and its 
backbone (see Fig. [TTl). 

We have studied RWs both on the percolation cluster and its backbone, performing 
10^ realizations on each cluster and average over 1000 clusters in each space dimensions 
d = 2,3,4. Estimates of scaling exponents u^ and I'^y^, describing scaling of walks 
on percolation cluster and backbone, respectively, are obtained by linear least-square 
fitting and given in Table [5l One can see, that the inequality z/^^ < z/^^ holds, and 
that the quantitative difference between these two values grows with increasing the 
space dimension d. On the other hand, both values are smaller than the corresponding 
exponent z/rw = 1/2, governing scaling of random walks on the pure lattice: presence of 
disorder slows down the diffusion process. The reason for this subdiffusive behavior is 
intuitively clear: due to the presence of defects, the randomly walking particle returns 
back to already visited sites more often, thus its walking distance is shorter than on the 
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Figure 12. Disorder averaged distribution function rP{r, N) vs the scaling variable 
r/{r) in rf = 2 dimensions, left: percolation cluster, right: backbone of percolation 
cluster. Lattice size L — 400, number of RW steps A^=180 (squares), 7V=160 
(triangles), 7V=140 (crosses). 
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Figure 13. Disorder averaged distribution function rP{r, N) vs the scaling variable 
r/(r) in d = 3 dimensions, left: percolation cluster, right: backbone of percolation 
cluster. Lattice size L = 200, number of RW steps A^^lOO (squares), N=90 (triangles), 
iV=80 (crosses). 
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Figure 14. Disorder averaged distribution function rP{r, N) vs the scaling variable 
r/(r) in d = 4 dimensions, left: percolation cluster, right: backbone of percolation 
cluster. Lattice size L = 50, number of RW steps N=60 (squares), N=50 (triangles), 
N—40 (crosses). 



pure lattice. This has also been observed in recent studies of less disordered deterministic 
fractals such as two-dimensional Sierpinski carpet composites [76] . 

The disorder averaged distribution function, defined in Eq. ([8]), and rewritten in 
terms of the scaled variables r/{r) as: 

rP(^^ ~ /(r/^) ~ fir/N'""') (11) 

is shown in Figs. [121 [IS and [14] for d = 2,3, 4, both for the cases of percolation cluster 
and backbone. When plotted against the scaling variable r/N'^^", the data are indeed 
found to nicely collapse onto a single curve, using our values for the exponent v^, i^rw 
reported in Table [3 
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5. Conclusions 

We studied the scaling behavior of simple random walks and self-avoiding walks on 
disordered lattices. Both models are of great interest: RWs provide a good description 
of diffusion processes, SAWs are successful in describing the universal properties of long 
flexible polymer macromolecules in a good solvent. 

We consider the case, when concentration p of lattice sites allowed for walking 
equals the critical concentration pc and the incipient percolation cluster can be found 
in the system. Studying properties of percolative lattices, one encounters two possible 
statistical averages: in the first, one considers only percolation clusters with linear size 
much larger than the typical length of the physical phenomenon under discussion, the 
other statistical ensemble includes all the clusters, which can be found on a percolative 
lattice. In our study, we considered only the first case, being interested in random 
and self-avoiding walks on a percolation cluster, which has a fractal structure. In this 
regime, the critical behavior of both models belongs to a new universality class: scaling 
law exponents change with the dimension of the (fractal) lattice on which the walk 
resides. 

We performed numerical simulations of random and self-avoiding walks on 
percolation clusters and the backbone of percolation clusters on lattice sizes L = 
400, 200, 50 in space dimensions d = 2,3, 4, respectively. Our results bring about the 
estimates for critical exponents, governing the scaling behavior of the models. We 
found that the statistics of SAWs is governed by the same scaling exponent both on 



Table 5. Exponents governing the scaling law of the end-to-end distance for RWs 
on percolation clusters and the backbone of percolation clusters. RS: real-space 
renormalization-group, EE: exact enumerations, MC: Monte Carlo simulations. 



^'rW\c^ 


2 


3 


4 


RS t20j 


0.356 


0.285 




EE [25 


0.349±0.002 






m 




0.266±0.01 




analytic [23j 


0.352 


0.268 




MC [25 


0.352 ±0.006 


0.271 ±0.004 




m 


0.352 ±0.006 


0.271 ±0.004 




m 


0.392 ±0.007 


0.282 ±0.003 




m 


0.348 ± 0.009 


0.274 ±0.007 




m 






0.222 ±0.007 


our results 


0.353 ±0.003 


0.273 ±0.003 


0.231 ±0.003 


^RW\'^ 


2 


3 


4 


analytic [5S] 


0.371 ±0.001 






EE[30] 


0.372 ±0.005 






MC [S] 


0.370 ±0.003 






our results 


0.372 ±0.002 


0.306 ±0.002 


0.262 ±0.002 
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a percolation cluster and its backbone: since the walks, which visit the dead ends 
are eventually terminated after a certain number of steps, the walks that survive in 
the limit A^ ^ oo on a percolation cluster are those confined within its backbone. For 
simple random walks, however, the picture is different: they cannot be trapped in "dead 
ends" . However, visiting these parts of a cluster will lead to some "slowing down" of the 
diffusion process in comparison with the behavior on the backbone where all the dead 
ends are removed. We found that the inequality u^ < v^ holds, and the quantitative 
difference between these two values grows with increasing space dimension d. 

The presence of disorder leads to a stretching of the trajectory of self-avoiding 
walks: the value of t'sAw i^ larger than the corresponding exponent on the pure lattice. 
However, the exponent v^^-, governing scaling of random walks on percolative lattices is 
smaller than that on a pure lattice: presence of disorder slows down the diffusion process. 
This can be explained as follows: due to the presence of defects, the randomly walking 
particle returns back to already visited sites more often, thus its walking distance is 
shorter than on the pure lattice. 
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7. Appendix 

To estimate the critical exponents t'sAw Tsaw linear least-square fits with varying lower 
cutoff for the number of steps A^min are used in order to detect possible corrections to 
scaling. For estimating z^saw ^® ^^^ hnear fits for the averaged end-to-end distance 
ln((r(A^))) = ln(y4) + z/g^-^^lnA^, and for 7gAw '^^ employ Eq. (ITOl) . Since this is an 
important aspect for assessing the quality of our final exponent estimates discussed 
in the main text, we have compiled in this Appendix these more detailed results in 
Tables IMTTl The x^ value (sum of squares of normalized deviation from the regression 
line) divided by the number of degrees of freedom, DF, given in the last rows, serves as 
a test of the goodness of fit. 
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Table 6. Results of linear fitting of obtained results for (r) for SAWs in d=2 
dimensions on the backbone of percolation clusters, L—4QQ. 



A^min 


'^SAW 


A 


XVDF 


11 


0.790 ± 0.005 


0.829 ± 0.003 


2.396 


16 


0.786 ± 0.005 


0.841 ± 0.005 


1.910 


21 


0.782 ± 0.004 


0.847 ± 0.005 


1.479 


26 


0.783 ± 0.003 


0.842 ± 0.006 


1.262 


31 


0.782 ± 0.003 


0.840 ± 0.007 


0.839 



Table 7. Same as Table [6] for d=3, L=200. 



iVmin 


'^SAW 


A 


X'/DF 


11 


0.668 ± 0.003 


0.935 ± 0.004 


2.269 


16 


0.669 ± 0.003 


0.930 ± 0.004 


2.054 


21 


0.669 ± 0.003 


0.924 ± 0.004 


1.345 


26 


0.667 ± 0.002 


0.930 ± 0.006 


0.743 


31 


0.668 ± 0.002 


0.934 ± 0.008 


0.844 



Table 8. Same as Table H for d=4, L=50. 



A^min 


'^SAW 


A 


X'/DF 


8 


0.588 ± 0.002 


1.025 ± 0.004 


2.615 


10 


0.587 ± 0.002 


1.023 ± 0.006 


1.777 


12 


0.586 ± 0.003 


1.021 ± 0.01 


0.978 


14 


0.586 ± 0.003 


1.031 ± 0.01 


0.767 



Table 9. Results of linear fitting of obtained results for Cn for SAWs in d—2 
dimensions on the backbone of percolation clusters, L=400. 



7ilw A xVDF 



16 


1.341 ±0.005 


1.219 ±0.003 


3.135 


21 


1.349 ±0.005 


1.189 ±0.003 


2.682 


26 


1.351 ±0.007 


1.168 ±0.002 


1.913 


31 


1.352 ±0.008 


1.172 ±0.002 


1.621 


36 


1.350 ±0.008 


1.163 ±0.001 


0.704 
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Table 10. Same as Table H for d=3, L=200. 



11 


1.265 ± 0.004 


1.82 ± 0.003 


2.767 


16 


1.268 ± 0.005 


1.192 ± 0.003 


2.135 


21 


1.270 ± 0.006 


1.184 ± 0.003 


1.968 


26 


1.267 ±0.008 


1.176 ±0.002 


1.513 


31 


1.269 ±0.008 


1.172 ± 0.002 ± 


0.762 



Table 11. Same as Table |9] for d=4, L=50. 

8 1.251 ± 0.005 1.77 ± 0.003 1.767 
10 1.252 ± 0.007 1.182 ± 0.003 1.135 
12 1.250 ± 0.008 1.184 ± 0.003 0.968 
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